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SUMMARY 


Flutter dynamic pressure and flutter frequency were obtained for a 
rectangular unstressed isotropic panel with all edges clamped for the follow- 
ing sets of parameters: (a) Panel length-to-width ratio of 0.5^ Mach number 

range of I .05 to 1.40, and ratio of boundary-layer thickness to panel length 
of 0.032 to O.lllj (b) panel length-to-width ratio of 2.0, Mach number of 1.20, 
and ratio of boundary-layer thickness to panel length of O.OI 6 to O.O 56 . For 
the configurations tested, it is shown that the turbulent boundary layer has 
a large stabilizing influence on the flutter of flat panels . The effect on 
flutter dynamic pressure is maximum near M = 1.20 and decreases rapidly with 
increasing Mach number. The effect on flutter frequency is minimum near 
M = 1.20 and increases with increasing Mach number. 


INTRODUCTION 


There continues to be substantial disagreement between theoretical and 
experimental panel flutter results in the low supersonic speed range ( refs . 1 
and 2 ). It has been suggested (ref. 2), and some experiments have indicated 
(ref. 3 ), that part of the disagreement may be due to the presence of the tur- 
bulent boundary layer. The influence of the boundary layer on panel flutter 
has particular significance for large launch vehicles and large supersonic 
aircraft because the boundary -layer thickness on substantial portions of such 
vehicles may not be small relative to panel dimensions. Prior to the present 
investigation, there has been no systematic experimental evaluation of the 
effect of the boundary layer on panel flutter. 

In recent years a number of experimental studies (refs. 1-5) of the 
flutter of flat isotropic panels in low supersonic flows have been published. 
Experimental information on boundary-layer effects is contained in refer- 
ences 3 and 5 . Reference 3 shows that a small variation in boundary -layer 
thickness produces a significant change in flutter dynamic pressure; whereas 
reference 5 shows that a large change in boundary-layer thickness causes 
essentially no change in the flutter phenomena. However, in reference 5^ ^ 
number of experimental difficulties such as large temperature gradients in 
the panel and nonuniform static pressure distributions were noted which raise 
questions about the reliability of the results , 



All the above data show considerable scatter and substantial disagreement 
with theory. These disparities have been attributed to various sources in 
addition to the viscous boundary layer, among which are static pressure differ- 
ential, thermal stresses, built-in stresses, cavity effects, inaccurate or 
unknown panel edge conditions, nonuniform or unknoTO flow conditions, unknown 
structural damping, wind-tunnel-wall effects and imprecise methods of defining 
the flutter boundary. 

A number of theoretical analyses (refs. 6 - 9 ) of the effects of the 
boundary layer on panel flutter have been published. All have concluded that 
the boundary layer can significantly affect panel flutter, particularly in the 
low supersonic speed range; but unfortunately all of them suffer from various 
inadequacies, such as being restricted to two-dimensional panels, being 
restricted to infinite chord panels, being based on two-dimensional aerody- 
namics, or oversimplification of the boundary -layer model, and thus do not 
allow reasonable quantitative comparisons with experiment. 

The purpose of the present study was to eliminate as many as possible of 
the usual sources of experimental inaccuracy, and in particular, to obtain 
quantitative information concerning the influence of the turbulent boundary 
layer on panel flutter in the low supersonic speed range. 

The flutter boundary was determined for an unstressed rectangular panel 
with clamped edge conditions for ratios of length to width of 0-5 over the 
Mach number range of 1.05 "to lAO and 2.0 at M = 1.20. The geometric 
boundary- layer thickness was varied from approximately O .25 to 1.0 inch. 
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P^a 

[ig cavity mass ratio parameter, 

V Poisson's ratio 

p density^ slugs/ft^ 

T panel thickness -to-length ratio, ^ 

Subscripts 

c cavity 

F flutter 

s structure (panel) 

00 free stream 


APPARATUS 


Panel Configuration 

Dimensional and mounting details of the panel and frame are shown in 
figure 1. The 9'00-inch by l8. 00-inch by 0 .0401-inch-thick panel was machined 
from a single piece of AZ31B-H24 magnesium alloy, 11.00 inch by 20.00 inch by 
0-100 inch. The active part of the panel was milled to final dimensions in a 
series of machining passes and was stress relieved after each pass to elimi- 
nate internal stresses and to maintain flatness. This construction produced 
a panel integral with a surrounding frame which was more than 15 times as 
stiff as the panel. This frame was attached to a secondary magnesium frame, 
1.0 inch wide by O. 5 O inch thick, with counterbored button head screws which 
produced a clamping force without introducing midplane stresses. The result- 
ing combined frame is more than 3300 times as stiff as the panel. This con- 
struction was used to closely simulate clamped edge conditions while 
maintaining a stress free panel - 

The panel frame was attached to the test fixture by one high stiffness 
and two low stiffness flexures as shown in figure 1. This three point mount- 
ing was designed to minimize the introduction of external stresses into the 
panel frame. 
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Wind Tunnel 


The experiments were conducted in the Ames 2- by 2-foot transonic wind 
tunnel, which is of the continuous flow type with porous test section walls 
surrounded by a plenum chamber (see ref. lO). Mach number is continuously 
variable over the range from 0.60 to 1.4o, and dynamic pressure is 
continuously variable over the range from 200 to I 5 OO psf . 


Variable Boundary-Layer Test Fixture 

The fixture for varying the boundary- layer thickness consisted of a 
splitter plate with a sharp leading edge installed in one side wall of the 
wind tunnel. This splitter plate could be positioned by remote control to be 
flush with the wall or to project into the airstream up to 1.0 inch. Dimen- 
sions of the test fixture are shown in figure 2, and front and rear view 
photographs are shown in figure 3* When the splitter plate is flush with the 
tunnel wall, the boundary- layer thickness in the test region is a maximum. 
When the splitter plate is moved into the airstream, a portion of the tunnel 
wall boundary layer is directed into the plenum chamber surrounding the test 
section and a thinner boundary layer is established on the splitter plate. 

The boundary- layer thickness in the panel test region is variable by this 
means between the approximate limits of O .25 and 1.00 inch. Four retractable 
boundary-layer total -pressure probes were employed to measure boundaiy-layer 
thickness at the points shown in figure 2. The variable boundary-layer test 
fixture containing the test panel is shown installed in the side wall of the 
2- by 2-foot transonic wind tunnel in figure 4. 


Other Equipment 

The panel was backed by a sealed rectangular cavity with the same 
dimensions as the panel and with an effective depth of 7*5 inch (where effec- 
tive depth is defined as the actual cavity volume minus volume occupied by 
instrumentation all divided by panel area). The panel frame was sealed to the 
cavity with adhesive plastic tape to minimize the possibility of introducing 
stress into the panel frame. 

The differential pressure between the cavity and the free stream was 
controlled manually to ± 0.25 psf by a system similar to that described in 
reference 1. 

The amplitude of flutter could be controlled by the pressure differential 
between the cavity and free stream. As an additional safety feature, mechani- 
cal "flutter stoppers" were installed in the panel cavity. The "flutter stop- 
pers" consisted of double-acting pneumatic cylinders topped with sponge rubber 
pads, which, when the cylinders were activated, pressed the pads against the 
panel surface and stopped the flutter motion mechanically. 
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PMEL NATURAL FREQUENCIES 


To determine the degree to which the panel conformed to the theoretical 
model of a homogeneous^ unstressed panel with clamped edge conditions ^ the 
natural frequencies and associated node lines were determined experimentally. 
Because the natural frequencies are influenced by many quantities, among which 
are the stiffness and density of the panel material, the density of the air 
surrounding the panel, the depth of the cavity behind the panel, and stresses 
such as those produced by a temperature difference between the panel and frame, 
the effect of each of these variables was determined. 


Panel Stiffness 

Panel stiffness was determined by measuring the deflection in the center 
of the panel with a noncontacting capacitive type displacement transducer 
while applying a uniformly distributed load with differential air pressure and 
equations presented in reference 11. The modulus of elasticity of 
5*6X10 psi was computed on the basis of average measured panel thickness of 
0.0^01 inch and an assumed Poisson’s ratio of 0.35* Panel thickness was deter- 
mined with an ultrasonic thickness measuring instrument and material density 
was determined by weighing a sample in air and in water. 


Air Density Effects 

The effect of air density on natural frequencies, with the panel backed 
by a 7 • 5 ”inch-deep cavity essentially identical to the one used during the 
wind-tunnel test, is shown in figure 5 ty the change of panel natural frequen- 
cies with air pressure. These natural frequencies can be extrapolated to zero 
pressure where there will be no air density or cavity effects . In the follow- 
ing table these frequencies are compared with theoretical in vacuo values 
determined by a 196 -mode (l4 modes in each direction) Galerkin solution utiliz- 
ing clamped beam modes (ref. 12 ). The measured node lines associated with 
each natural frequency are also shown. This comparison indicates that the 
panel confoms within 1 percent to the theoretical model for the first two 
natural frequencies. It will be noted that the mode shapes for the fourth and 
fifth measured natural frequencies are in reverse order from theory. This 
switch in mode shape sequence can easily be accounted for by an approximately 
5 percent difference in stiffness between the panel axes . This difference in 
stiffness could result from the use of rolled material for fabricating the 
panel as well as from the panel support system which is more rigid along the 
short panel axes . Because the primary interest in this investigation is in 
the first two panel natural frequencies, it can be said that the panel 
conforms very well to the above theoretical model. 
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Panel Natural Frequencies and Locations of Node Lines 


Location of 
node lines 


1 

1 
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1 

1 

1 1 

1 1 

i 1 

« * 


s * ! 

1 1 1 

1 1 1 

« . 

Measured natural 
f requeue ies , Hz 

110 

i43 

212 

287 

298 

Computed natural 
f requeue ies , Hz 

109.7 

14-2.0 

199.8 

285.5 

282.6 

Deviation of 
measured natural 
frequencies from 
computed values, 
percent 

0.27 

i 

0.70 

6.11 
1 

0.53 

5.45 


Cavity Depth Effects 

The effect of cavity depth on the natural frequencies of the unstressed 
panel at atmospheric pressure is shown in figure 6, These data were obtained 
with an unvented variable depth cavity with all walls more than 1000 times as 
stiff as the panel. These results indicate that a 7 ^^-Inch-d-eep cavity at 
atmospheric pressure will Increase the first natural frequency approximately 
9.4 percent above the in vacuo natural frequency, but will change the higher 
natural frequencies less than 2 percent. Since all flutter data presented in 
this report were taken at static pressures between O.25 and O.65 atmospheres, 
the influence of the cavity on these data will be less than indicated above. 
Figure 6 also indicates that for a very shallow cavity, the natural frequency 
for the first mode is higher than for the second mode, which agrees with the 
information presented in reference 13* The natural frequencies presented in 
figures 5 and 6 were determined by the forced vibration technique of 
reference l4. 


Thermal Effects 

The theoretical effect on panel natural frequencies of a uniform 
temperature differential between the panel and frame is shown in figure 7 for 
compressive resultant stress. The temperature at which the first natural fre- 
quency goes to zero corresponds to buckling of the panel (ref. 15 ); in this 
case 3*3*^ F. This figure indicates that small temperature differentials can 
have a large effect on panel natural frequencies and, consequently, a large 
effect on the flutter boundary. 

In actual practice it is not possible to get a uniform temperature 
differential because of thermal conduction between the panel and frame; thus, 
an experiment was conducted to determine the actual effect of temperature 
differentials on natural frequencies. The pa^^i radiantly heated very 
rapidly and its natural frequencies were determined for various temperature 
differentials. In an attempt to conform as nearly as possible to theory, the 
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frame was insulated from the heat source; and heat was concentrated near the 
panel houndaries to compensate for conduction into the frame. Typical temper- 
ature distributions measured on the unheated side of the panel for various 
times from start of heating are shown in figure 8. The panel natural frequen- 
cies were determined by spectral analysis of the transient displacement signal 
resulting from periodically exciting the panel with a mechanical impulse dur- 
ing the heating cycle. The natural frequencies as a function of temperature 
differential between the frame and the center of the panel are shown in fig- 
ure 9 • Buckling was observed at a differential temperature of approximately 
12^ F which coincides with the differential temperature for minimum resonant 
frequencies. The difference between theory and experiment in buckling temper- 
ature is probably due to the inability to obtain a step change in temperature 
at the boundary because of thermal conduction and also to possible panel cui^a- 
ture due to temperature differentials through the panel caused by heating on 
only one side (ref. 15 )- 

The above experiment indicates that a ±0.5^ F temperature differential 
will change the first natural frequency less than 2 percent. Thus it is 
expected that a temperature differential of this magnitude would have a corre- 
spondingly small effect on the flutter boundary (ref. l6). Consequently; the 
panel-frame temperature differential was maintained at less than ±0.5*^ F during 
the wind-tunnel test; with most of the variation being less than ±0.2^ F. 


INSTRUMENTATION 


Panel displacement was measured with a noncontacting capacitance type 
transducer which produced an output voltage linearly proportional to the panel 
displacement. This instrument; which has a range of 0.100 inch; a flat fre- 
quency response of dc to 10;000 Hz, and linearity of ±1.0 percent; was posi- 
tioned 0.050 inch from the under surface of the panel at the location shown in 
figure 2. 

The steady surface static pressures and the total pressures from the 
boundary-layer probes were measured with conventional strain gage differential- 
pressure transducers connected to the orifices by short lengths of tubing and 
referenced to free-stream static pressure. 

The unsteady surface pressures were measured with flush diaphragm pressure 
transducers with piezoresistive sensing elements. These transducers have a 
diameter of O.I25 inch; a pressure range of ±1.0 psi differential; and a flat 
frequency response of 0 to 12;000 Hz, and were referenced to a local surface 
static pressure through a length of small diameter tubing to damp frequencies 
above 10 Hz. 

Panel and frame temperatures were measured with iron-constantan 
thermocouples in the center of the panel and at two points on the frame . The 
thermocouple in the center of the panel was constructed of wire O.OO5 inch in 
diameter and was the only object contacting the panel. 
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steady qiiantities such as temperatures and pressures were recorded 
digitally. Unsteady or slowly varying quantities such as fluctuating static 
pressure^ panel response, and pressure differential between the cavity and 
free stream were recorded by a frequency modulated magnetic tape recorder with 
a frequency response of dc to I 25 O Hz, 


TEST PROCEDURE 


Flow Conditions in Panel Test Region 

An important consideration in panel flutter testing is that the sirrface 
pressure distribution in the test region be uniform. A nonuniform static 
pressure distribution would cause local differential pressures between the 
cavity behind the panel and the free stream which would stress and distort the 
panel locally, thereby producing erroneous flutter data. Typical values of 
the surface static pressure coefficient in the test region for the present 
tests are shown in figure 10. The maximum variation is ±0.02, with most 
variations being less than ±0.01. 

An additional consideration in panel flutter testing is that fluctuating 
surface pressures such as those produced by boundary- layer turbulence, frec- 
stream turbulence, and wind-tunnel compressor drive harmonics not contain any 
strong sinusoidal components in the frequency range where flutter occurs. 

Such components could cause the panel to exhibit flutter-like response under 
low damping conditions prior to actml flutter. Representative values of the 
root-mean-square fluctuating pressure coefficient at M = 1.20 for various 
boundary-layer thicknesses at stations upstream and downstream of the panel 
test region are shown in figiire 11. A typical power spectrum of the fluctuat- 
ing surface pressure is shown in figure 12. This power spectrum was obtained 
with a constant bandwidth (5 Hz measured at the 3 points) spectrum analyzer 
with linear sweep from a data sample of approximately 20-seconds duration. 

The power spectrum indicates that there are no dominant frequencies in the 
frequency range of Interest of the present investigation. 

A typical variation of the geometric boxuidary- layer thickness with 
splitter plate position is shown in figure 13* Typical velocity profiles 
along the wall centerline for various splitter plate displacements are shown 
in figures l4(a) and l4(b) for stations upstream and downstream, respectively, 
of the panel test region. A l/7th power profile is shown in both figures for 
comparison . 

The effect of the wind-tunnel wall opposite the panel on the flutter 
boundary will be negligible from purely geometric considerations because 
acoustic disturbances originating on the panel are reflected downstream for 
all Mach number and configuration combinations tested- The effects of the 
wind-tunnel side walls are expected to be small because of geometric consider- 
ations as well as the large attenuating effects of the porous tunnel walls and 
the turbulent boundary layer on these walls as discussed in reference 17 • 
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Determination of the Flutter Boundary 

The criterion for a flutter boundary in a theoretical analysis is very 
simple; however, the experimental determination of this boundary has proved to 
be very difficult (ref. l). 

From a theoretical viewpoint the flutter boundary is defined as the locus 
of neutral stability. For a linear theory the stable side of a flutter bound- 
ary corresponds to decaying initial disturbances and the unstable side to 
imbounded an^lif ication of initial disturbances. For a nonlinear theory, the 
stable side of a flutter boundary again corresponds to decaying initial dis- 
turbances but the unstable side is characterized by a limit cycle oscillation 
for which the amplitude goes to zero as the boimdary is approached from the 
unstable side. 

The experimental determination of the flutter boundary is complicated by 
several factors among which are the nonlinear behavior of clamped plates for 
finite deflections producing in-plane tension, the presence of a random 
\insteady pressure field due to boundary layer and free-streara turbulence, 
slowly changing free-stream conditions, and variation in pressure differential 
between the cavity and free stream due to the changing free-streara conditions. 
The random unsteady pressure field produces forced random motion of the panel 
with largest amplitudes in the frequency spectrum in the vicinity of the nat- 
ural frequencies of the panel-airstream system. As the flutter boundary is 
approached from the stable side, the total damping of the system decreases at 
frequencies near the flutter frequency (ref. 1^); consequently, panel ampli- 
tude due to the fluctuating pressures will Increase. At or near the flutter 
boundary, the panel-airstream system acts like a narrow band-pass filter with 
the random pressure field the input and the panel response the output. At the 
flutter boundary, the limit cycle motion due to flutter is of zero amplitude 
and so undetectable. As the unstable region is penetrated, the limit cycle 
grows in amplitude and eventually dominates the random response. Now the out- 
put of a narrow band-pass filter with a random input is hard to distinguish 
from a sinusoid (ref. l8); thus it is difficult to separate the increased ran- 
dom response due to decreased damping from the true flutter response. 

The criteria used in past works to define the flutter boundary 
experimentally have either not been stated at all or only imprecisely . The 
problem is discussed thoroughly in reference 1 where results are presented in 
the form of a scatter band for each flutter point, i.e., an upper and lower 
bound for" each flutter point is stated. While this is perhaps the most satis- 
factory method’ used to date for presenting experimental results, it leaves the 
flutter boundary uncertain within the width of the scatter band, and, of 
course, there remains arbitrariness in selection of the upper and lower limits . 

In the present study, the flutter boundary is determined as follows ; At 
dynamic pressures well below that for flutter, the panel response to pressure 
disturbances from the turbulent boundary layer increases approximately lin- 
early with dynamic pressure, as indicated in figure 15- In this figure, panel 
amplitude in the mode which becomes the flutter mode is shown as a function of 
dynamic pressure. This "subflutter response" is used as a reference level to 
define the flutter point systematically. The straight lines labeled 2x, 3x, 
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etc., in figure 15 are multiples of the subflutter response. Their intersec- 
tions with the amplitude response curve define values of dynamic pressure 
necessary to produce panel amplitudes which are these multiples of the sub- 
flutter response. Figure l6 shows dynamic pressure as a function of the mul- 
tiples of the subflutter response for the same conditions as shown in fig- 
ure 15. As the multiple increases, the difference in dynamic pressure between 
successive multiples decreases; that is, the dynamic pressure appears to be 
converging to a limit. Furthermore it appears that the 6x multiple results in 
a dynamic pressure very close to this limiting dynamic pressure and is 
therefore arbitrarily defined herein as the flutter boundary. 

The following test procedure was adhered to throughout the program. Mach 
number and splitter plate position were held constant while dynamic pressure 
was increased in small Increments beginning at a value well below the flutter 
boundary. After each increase in dynamic pressure, the panel-frame tempera- 
ture difference was allowed to stabilize to ±0.5° F to minimize thermal stress 
effects. The indicated differential pressure between the cavity and free 
stream was then varied slowly from -50 psf to +50 psf while the panel response 
was recorded. It is assumed that the maximum response so measured corresponds 
to zero average differential pressure across the panel. Next the indicated 
differential pressure was set to the value at which maxlmiim response occurred. 
This response was recorded so that amplitude spectrums could be obtained with 
all test conditions fixed. The above procedure was repeated for each increase 
in dynamic pressure until the maximum response was a large ati^jlitude, essen- 
tially sinusoidal motion indicating that the flutter boundary had been pene- 
trated. Boundary-layer velocity profiles were obtained at dynamic pressures 
near the flutter boundary. The foivard probe was retracted while data were 
being recorded for the rear probes, and all probes were retracted when bound- 
ary layer data were not being recorded. 

Panel amplitude response at dynamic pressures below flutter was a random 
waveform in which most of the energy was concentrated near the natural fre- 
quencies of the panel-airs t ream system. The special data reduction system 
shown in block diagram form in figure 17 was developed to determine the ampli- 
tude response of the panel mode which becomes the flutter mode. This system 
automatically tracks and filters a dominant response peak as its frequency and 
amplitude vary with differential pressure. The system outputs voltages proper 
tional to the frequency it is tracking and to the amplitude at that frequency. 
These quantities are automatically plotted as functions of the differential 
pressure between the cavity and free stream. A representative plot produced 
by this system is shown in figure 18. Ebccept for the data in figure 19 ^ s-ll 
panel response and flutter frequency data presented in this report were 
obtained from such plots . 

As noted previously, panel response was also recorded at the fixed 
differential pressure between the cavity and free stream at which this 
response was maximum. An amplitude spectrum for such a fixed differential 
pressure is shown in figure 19. The spectrum was obtained with a constant 
bandwidth (5 Hz measured at the 3 dB points) spectrum analyzer with linear 
sweep from a data sample of approximately 20 seconds duration. This figure 
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corresponds to the same test conditions as those in figure l8 with which it 
agrees well in hoth frequency and amplitude. 


RESULTS AND DISCUSSION 


Influence of the Turbulent Boundary Layer on Flutter Dynamic Pressure 

The dynamic pressures for various multiples of subflutter response for a 
series of boundary- layer thicknesses at M = 1.20 are shown in figure 20 for 
a/b =0.5 and in figure 21 for a/b = 2.0. The value of the geometric boundary- 
layer thickness used is the thickness in the center of the panel estimated by 
assuming linear growth between the upstream and downstream probes. The curve 
for 6x subflutter response is considered the flutter boundary. The flutter 
dynamic pressure in figures 20 and 21 is approximately a linear function of 
the geometric boundary-layer thickness . 

In figures 20 and 21 the rate of change of flutter dynamic pressure with 
boundary- layer thickness is much greater for the a/b =0.5 case where the 
panel length is 9*0 inches than for the a/b =2.0 case where the panel length 
is l8.0 inches. If the boundary -layer thickness is normalized with respect to 
the panel length (b/a)^ as discussed in reference 7^ the rate of change of 
flutter dynamic pressure is much more nearly the same for the two different 
panel lengths as shown in figure 22^ thus indicating that an important param- 
eter in the study of boundary-layer effects is b/a^ However, it should be 
noted that a/b is not identical; therefore, the two curves would not be 
expected to be identical. 

The dynamic pressure for various multiples of panel subflutter response 
as a function of Mach number for a/b = 0-5 is shown in figure 23 for the maxi- 
mum geometric boundary-layer thickness (approximately 1.0 in.) and in fig- 
ure 24 for the minimum geometric boundary- layer thickness (approximately 
0.25 in.). In both cases the boundary -layer thickness varies slightly with 
Mach number. Both figures show the characteristic minimum (refs. 1, 2, 3) in 
flutter dynamic pressure in the vicinity of M = 1.20 and a rapid increase with 
increasing Mach number. 

The approximate linear relationship of flutter dynamic pressure with 
boundary- layer thickness indicated in figures 20 and 21 is used to show the 
change of flutter dynamic pressure with boundary-layer thickness (dq/db) for 
a/b = 0.5 as a function of Mach number in figure 25. This figure indicates 
that the maximum boundary- layer effect occurs near M = 1.20 and decreases 
rapidly with increasing Mach number. 

By use of a linear relationship of flutter dynamic pressure with boundary 
layer thickness at all Mach numbers, these data can be reasonably extrapolated 
to zero boundary- layer thickness to obtain a flutter boundary for comparison 
with results predicted by potential flow theory. This linear relationship can 
also be used to derive the flutter dynamic pressures at various other constant 
boundary-layer thicknesses (for a/b = 0.5) shown in figure 26. 
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Influence of the Turbulent Boundary Layer on Flutter Frequency 

The flutter frequencies for various boundary -layer thicknesses at a Mach 
number of 1.20 are shown for a/b = 0.5 snd 2.0 in figure 27. Here a linear 
relationship with geometric boundary-layer thickness is again noted. 

The rate of change of flutter frequency with boundary- layer thickness 
(df/dB) at M = 1.20 is much greater for a/b = 0.5 (panel length of 9.0 in.) 
than for a/b = 2.0 (panel length of l8.0 in.) as indicated in figure 27- Nor- 
malizing the boundary- layer thickness with respect to the panel length (B/a) 
makes the rate of change of flutter frequency with this quantity approximately 
the same for both values of a/b as shown in figure 28. 

Using a linear relationship of flutter frequency to boundary- layer 
thickness for all Mach numbers, figure 29 is derived showing flutter frequency 
as a function of Mach number for various constant geometric boundary-layer 
thicknesses for a/b =0.5. There is some indication from the figure (for 
B = 0.95 inch at least) that the flutter frequency approaches zero as Mach 
nimber decreases toward 1.0. This is in qualitative agreement with the 
theoretical results presented in reference 19- 

The rapidly increasing flutter frequency in the vicinity of M = 1.40 for 
thin boundary layers is probably associated with the change in the nature of 
the flutter from single degree of freedom to two mode coupled type as pre- 
dicted by theory in this Mach number range. Based on results of reference I9, 
it is expected that for a thin boundary layer, the flutter frequency versus 
Mach number curve will again level out near M = I.60, after which it will 
change slowly. For a thick boundary layer the rapid change in frequency would 
be expected to occur at a higher Mach number. 

The rate of change of flutter frequency with boundary-layer thickness 
(df/dB) as a function of Mach number for a/b = 0.5 is shown in figure 30. The 
Increasing effect of boundary-layer thickness on flutter frequency in the 
range M = 1.20 to 1.4o is probably associated with the rapid change of flutter 
frequency with Mach number in this region while the Increasing effect of 
boimdary-layer thickness below M = 1.20 is associated with the rapid decrease 
of frequency toward zero as Mach number decreases to M = 1.0. It is interest- 
ing to note from figures 24 and 30, that the boundary layer has minimum effect 
on flutter frequency and maximtun effect on flutter dynamic pressure at approx- 
imately the same Mach number of 1.20. From these results, it is clearly evi- 
dent that Ignoring boundary- layer effects may lead to considerable difference 
between theory and experiment for flutter frequency and/or flutter dynamic 
pressure at low to moderate supersonic Mach numbers . 


Panel Flutter Parameters 

A natural set of dimensionless parameters for the flutter of flat 
isotropic panels exposed to linearized three-dimensional unsteady potential 
flow (ref. 20) are the mass ratio parameter, p = p^a/p^h, the stiffness param- 
eter, k =^J T^E/coo^Pg(l - v^), the panel length-to-width ratio, a/b, the Mach 
number M, and the reduced frequency k . Flutter boundaries may be plotted in 
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the M-k plane for various reduced frequencies for each Mach number and 
length-to-width ratio. In figure 31 the mass ratio parameter is plotted as a 
function of boundary- layer thickness for the various Mach nimibers . Since 
these data fall approximately on a straight line, they can reasonably be 
extrapolated to zero boundaiy- layer thickness to obtain a value of \i for 
comparison vith results predicted by potential flow theory. Again the influ- 
ence of the boundary-layer thickness is clearly evident. The reduced fre- 
quency as a function of Mach number for various constant geometric boundary- 
layer thicknesses is presented in figure 32. Computed values of and k and 
other pertinent quantities^ such as boundary- layer thickness and dynamic 
pressure for flutter^ are presented for each flutter point in table I- 


CONCLUSIONS 


1. The turbulent boundary layer has a large stabilizing influence on the 
flutter of flat isotropic panels at low supersonic Mach numbers . 

2 . The effect of the turbulent boundary layer on flutter dynamic 
pressure is maximum near M = 1.2 and decreases rapidly with increasing Mach 
number for the configurations tested. 

3 . The effect of the turbulent boundary layer on flutter frequency is a 
minimum near M = 1.2 and increases with increasing Mach number over the range 
M = 1.2 to 1.4 for the configurations tested. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 9^035, Dec. 5, I967 
126-14-02-01-00-21 
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TABLE I.- TABULATED PANEL FLUTTER DATA 
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Figure !•- Panel dimensions and construction details. 
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Figvire 3«- Variable boundary -layer test fixture* 
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(b) Three -foiirths rear viev of plenum chamber side 
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transonic vind tunnel. 




Figure 5»“ Effect of static pressure on the natural freq^uencies of the test panel (cavity depth 




Cavity depth, in. 

Figure 6.- Effect of cavity depth on panel natural frequencies at atmos 

pheric pressure . 
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Figure 7**" Theoretical panel natural frequencies as function of differential temperature (no 
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Figure 9-- Experimental panel natiiral frequencies as a function of differential temperature (no 
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Figure 10.- Typical wall centerline static -pres sure distributions in region occupi' 
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Figure I5.- Typical variation of maximum panel response with dynamic pressure 





Multiples of subflutter response 

Figure 16.- Typical convergence of multiples of subflutter response with 

dynamic pressure. 





Fig\H*e 17.- Block diagram of system used to reduce panel response data. 
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6, in. 

Figure 20.- Dynamic pressure for multiples of subflutter response as a func 
tion of boundary-layer thickness for length-to-width ratio of 0.5* 








Figure 22.- Flutter dynamic pressure as a function of the ratio of bound 

ary- layer thickness to panel length. 
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Figure 23 .- Dynamic pressure for multiples of subflutter response as a 
function of Mach number for maximuin boiuidary-layer thickness (approx- 
imately 1.0 in.). 
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Figure 24.- Dynamic pressure for multiples of subflutter response as a 
function of Mach number for minimum boundary -layer thickness (approx 
imately 0.25 in.). 




Figure 25*- Rate of change of flutter dynamic pressure with houndary-layer 
thickness as a function of Mach number- 
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Figure 26.- Flutter dynamic pressure as a function of Mach number derived 
from linear relationship of flutter dynamic pressure with boundary- 
layer thickness. 
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Figure 27* - Flutter frequency as a function of boundary -layer thickness 
for length-to-width ratios of O .5 and 2.0, 














Figure 31 •- Mass ratio parameter as a function of boundary -layer thickness 
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